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Introduction 


Nonlinear supersymmetry is a natural generalization of the usual linear supersymmetry HI- 
It is realized variously in such different systems as the parabosonic 0 and parafermionic || 
oscillator models, and the P, T-invariant models of planar fermions 0 and Chern-Simons 
fields [pj]. It is also the symmetry of the fermion-monopole system 0.0- The algebraic 
structure of the nonlinear supersymmetry resembles the structure of the finite W-algebras 
[T)j for which the commutator of generating elements is proportional to a finite order poly¬ 
nomial in them. In the simplest case the nonlinear supersymmetry is characterized by the 
superalgebra of the form 

[Q ± ,H\ = 0, (Q ± ) 2 — 0. {Q + ,Q-} = Pn(H), (l) 


where P n (.) is a polynomial of the n-tli degree. The nonlinear supersymmetry with such a 


superalgebra was investigated for the first time by Andrianov, Ioffe and Spiridonov |1C 
The pseudoclassical construction underlies the supersymmetric quantum mechanics of 
Witten E.0 corresponding to the linear (n = 1) case of the superalgebra (|1|). Though the 
nonlinear supersymmetry can also be realized classically, there is an essential difference from 
the linear case: the attempt to quantize the nonlinear supersymmetry immediately faces 

It was shown 11 that the universal algebraic 
in the form of the Dolan-Grady relations 
T3| underlies the so called holomorphic nonlinear supersymmetry [Tip. This structure allows 
ones to find a broad class of anomaly-free quantum mechanical systems related to the exactly 
and quasi-exactly solvable systems | [i~4| . [T~o| . 0 0, 0, [T9| , and gives a nontrivial centrally 
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the problem of the quantum anomaly 
structure with associated “integrability conditions 


extended generalization of the superalgebra O |J 


In this talk, after a short discussion of the intimate relation between the generalized 
statistics and supersymmetry ||, we shall review the recent results on the nonlinear super- 
symmetry obtained in the context of the quantum anomaly problem and of the universal 
algebraic construction associated with the holomorphic nonlinear supersymmetry 0 0 . 


Nonlinear supersymmetry in purely parabosonic 
systems 


Some time ago it was shown that the linear supersymmetry can be realized without fermions 
B 0 0- The nonlinear supersymmetry admits a similar realization revealing the close 
relationship between the generalized statistics and supersymmetry || . The relationship can 
be observed in the following way. Let us consider a single-mode paraboson system defined 
by the relations 

[{a + , a - }, a ± ] = ±a ± , a _ a + |0) = p|0), a“|0) = 0, 


where p G N is the order of a paraboson [23|. Then the direct calculation shows that the 
pure parabosonic system of the even order p = 2 [k + 1), k G Z + , with the Hamiltonian 
of the simplest quadratic form H = a + a~ reveals a spectrum typical for the nonlinear 
supersymmetry: all its states are paired in doublets except the k + 1 singlet states |2 1) oc 
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(a + ) 2, |0), l — 0,... ,k. In correspondence with this property, the system has two integrals 
of motion 

Q + — (a + ) 2fe+1 sin 2 ^{a + , a - }, = (a _ ) 2fc+1 cos 2 ^{a + , a - }, (2) 

which together with the Hamiltonian form the nonlinear superalgebra of the order n = 
2k + 1 with P 2k+ i(H) = H • nLi(^ 2 — 4m 2 ). This simplest system reflects the peculiar 
feature of the parabosonic realization of supersymmetry: the supercharges are realized in 
the form of the infinite series in a , and the role of the grading operator is played here by 
R = (—1)^ = costtN, where N = 4{a + , a~} — \p is the parabosonic number operator. 

It is known that the deformed Heisenberg algebra with reflection 

[a - , a + ] = 1 + uR, {i?, a ± } = 0, R 2 — 1, (3) 


underlies the parabosons [^, [25j. This algebra possesses unitary infinite-dimensional repre¬ 
sentations for v > —1, and at the integer values of the deformation parameter, v = p — 1, 
p € N, is directly related to parabosons of order p |23], ^4], On the other hand, at 
v = —(2 p + 1) the .R-deformed Heisenberg algebra has finite-dimensional representations 
corresponding to the deformed parafermions of order 2 p [^5|]. In the coordinate represen¬ 
tation the operator R is the parity operator and the operators ar can be realized in the 
form a ± = -m(x =F iD v 


t ^v) with D v = — — t^R). In the context of the Calogero-like mod¬ 

els, the operator D v is known as the Yang-Dunkl operator |26|, ^7|, where R is treated as 
the exchange operator. In the coordinate representation the Hamiltonian H = a + a~ and 
supercharges (|j) read as [3] 


H = 


d 2 


dx' 2 


+ x z + 


4a: 2 


Q + = W“)' = 


23 (fc+j) 


d v 

ia +x+ 2i 


- 1 HR 


(1 - R) 


2fc+l 


( 4 ) 


( 5 ) 


with v = 2k T 1. The system given by the Hamiltonian ([J) can be treated as a 2-pa.rticle 
Calogero-like model with exchange interaction, where x has a sense of a relative coordinate 


and R has to be understood as the exchange operator . Therefore, at odd values of the 

parameter u, the class of Calogero-like systems (||) possesses a hidden supersymmetry, which 
at v = 1 is the linear (n = 1) supersymmetry in the unbroken phase, whereas at v — 2k + 1, 
the supersymmetry is characterized by the supercharges being differential operators of order 
2k + 1 satisfying the nonlinear superalgebra (|). Recently the realization of the nonlinear 
supersymmetry was extended within the standard approach with fermion degrees of freedom 
to the case of multi-particle Calogero and related models |l0 . 


Classical supersymmetry 

Let us turn now to the classical formulation of the supersymmetry ([lj). For the purpose, we 
consider a non-relativistic particle in one dimension described by the Lagrangian 

L — -x 2 V(x) — L(x)N + i9 + 9~, (6) 


3 





















where 9 ± are the Grassman variables, (6* + )* = 9~, N = 9 + 9~, and V(x) and L(x) are two 
real functions. The nontrivial Poisson-Dirac brackets for the system are {x, p }* = 1 and 
{9 + , 9~} if = —i, and the Hamiltonian is 

H = l -p 2 + V(x) + L(x)N. (7) 

The Hamiltonian H and the nilpotent quantity N are the even integrals of motion for any 
choice of the functions V(x) and L(x), and one can put the question: when the system 
(||) has also local in time odd integrals of motion of the form Q ± = B T (x, p)9 ± , where 
(. B + )* = B~1 It is obvious that such odd integrals can exist only for a special choice of 
the functions V(x) and L(x). Restricting ourselves to the physically interesting class of the 
systems given by the potential V(x) bounded from below, we can generally represent it in 
terms of a superpotential W{x)\ V(x ) = \W 2 (x) + v, v G M. Then all the supersymmetric 
systems are separated into the three classes defined by the behaviour of the superpotential 
and the results can be summarized as follows [p~lJ1 . 

i) When the physical domain given by z = W ( x)+ip includes the origin z = 0 (a < W{x) < b, 
a < 0, b > 0), the corresponding supersymmetric system is characterized by the Hamiltonian 
and the supercharges of the form 

n 2 1 

H= 2 + 2 W ' 2 ^ + V + W '( x ) + W ( X ) M ( W2 ( X ))\ N , ( 8 ) 

Q+ = (( Q~y = z n e iSZM{p*-y*+W\x))dy e +, n £ 


where M(W 2 ) is an arbitrary regular function, |M(0)| < oo. The appearance of the inte¬ 
ger parameter illustrates in this case the known classical “quantization phenomenon” pi 


The appropriate canonical transformation reduces the system with these Hamiltonian and 
supercharges to the form of the supersymmetric system with the holomorphic supercharges 

0 : 


H = ip 2 + 1 W 2 (x ) + v + nW\x)9 + 9~, Q + = (Q~)* = z n 9 + , 
The integrals (|S|) obey the classical nonlinear superalgebra: 

{Q-, Q + }* = -iH n , {Q ± ,H}^ 0. 


n G Z + . 


(9) 


( 10 ) 


The presence of the integer number n in the Hamiltonian means that the instant frequencies 
of the oscillator-like odd, 9 ± , and even, z, z, variables are commensurable. Only in this case 
the regular odd integrals of motion can be constructed, and the factor z n in the supercharges 
Q ± corresponds to the n-fold conformal mapping of the complex plane (or the strip a < 
Rez < b) on itself (or on the corresponding region in C). 

ii) The physical domain is defined by the condition Re z > 0 (or Re z < 0) and also includes 
the origin of the complex plane. But unlike the previous case, there are no closed contours 
around z — 0. In this case the most general form of the Hamiltonian and the supercharge is 

H = Y + ^W 2 (x)+ v + w '(x)[a + R (W(x))\ N , Q+ = z a e if Z o^ cosA ) rf V, (11) 
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where a E M, and we assume that the function R{W) is analytical at W = 0 and -R(O) = 0. 


By the canonical transformation [|Il]] , the Hamiltonian and the supercharges can be reduced 
to 

H = ^p 2 + ^W 2 {x)+v + aW\x)6 + 6-, Q + = (Q~)* = z a 9 + , a E M+. 

iii) The physical domain is defined by the condition Rez > 0 (or Rez < 0), i.e. the origin of 
the complex plane is not included. Though in this case the Hamiltonian and the supercharges 
have a general form 


H = ^ + ^W 2 + v + W'(x)<l)(W{x))N, Q + = (Q~y = f(H)e if ^^ pcosX)dX 9 + , (12) 


where 0 is some function, the appropriate canonical transformation reduces it to [11 


H=\p 2 + \w 2 (x)+v, Q ± =e ± . 

This means that classically the supersymmetry of any system with bounded non-vanishing 
superpotential has a “hctive” nature. 

In what follows we will refer to the nonlinear supersymmetry generated by the holomor- 
phic supercharges with the Poisson bracket (anticommutator) being proportional to the n-th 
order polynomial in the Hamiltonian as to the holomorphic n-supersymmetry. 

Though the form of the Hamiltonians and ([12]) can be simplified by applying in 

every case the appropriate canonical transformation reducing the associated supercharges to 
the holomorphic or antiholomorphic form, the quantization breaks the equivalence between 
the corresponding classical systems (even in the linear case n = 1) [|ll|]. Moreover, alternative 


classical forms for the Hamiltonians and associated supercharges are important because of 
the quantum anomaly problem to be discussed below. Having in mind the importance 
of alternative classical formulations of the nonlinear supersymmetry from the viewpoint of 
subsequent quantization, one can look for the classical formulation characterized by the 
supercharges of the n-th degree polynomial form in p |1]]]. The problem of hireling such a 


formulation can be solved completely in the simplest case n = 2, for which the supercharges 
are given by 




while the Hamiltonian is 


( ± ip + W(x)) 2 + 


W 2 (x) 


e ± , 


c E 


J7 = i 

2 


p 2 + W 2 (x) - 


W 2 (x) 


+ 2W\x)N + v. 


(13) 


(14) 


Note that the Hamiltonian G) has the Calogero-like form: at W(x) = x its projection to 
the unit of Grassman algebra takes the form of the Hamiltonian of the two-particle Calogero 
system. Depending on the value of the parameter c, classically the Calogero-like n = 2 
supersymmetric system ([b|) is symplectomorphic to the holomorphic n-supersymmetry with 
n = 0 (c > 0), n = 1 (c < 0) or n = 2 (c = 0) |TT| . 
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Quantum anomaly and quasi-exactly solvable (QES) 
systems 


According to the results on the supersymmetry in pure parabosonic systems, a priori one 
can not exclude the situation characterized by the supercharges to be the nonlocal operators 
represented in the form of some infinite series in the operator jC Since such nonlocal 
supercharges have to anticommute for some function of the Hamiltonian being a usual local 
differential operator of the second order, they have to possess a very peculiar structure. We 
restrict ourselves by the discussion of the supersymmetric systems with the supercharges 
being the differential operators of order n. Classically this corresponds to the system (Bl) 


with the holomorphic supercharges or to the system (14). 

In the simplest case of the superoscillator possessing the nonlinear n-supersymmetry 
and characterized by the holomorphic supercharges of the form (|9|) with W{x) = x, the 
form of the classical superalgebra {Qt^Qn} = is changed for {Q^^Qn} = H(H — 
h)(H — 2 h)... (H — h(n — 1)). Moreover, it was pointed out in [[| that for W{x) ^ ax + b 
a global quantum anomaly arises in a generic case: the direct quantum analogues of the 
superoscillators and the Hamiltonian do not commute, [Q±,H n ] ^ 0. Therefore, we arrive 
at the problem of looking for the classes of superpotentials and corresponding quantization 
prescriptions leading to the anomaly-free quantum n-supersymmetric systems. 

Let us begin with the quantum supercharges in the holomorphic form corresponding to 
the classical n-supersymmetry, 


nn± 


q ± = (A^) n e 


where 


A ± = T h^- + W(x). 
ax 


(15) 


Choosing the quantum Hamiltonian in the form (i7j), from the requirement of conservation 
of the supercharges, [Q , H] = 0, one arrives at the supersymmetric quantum system given 
by the Hamiltonian [|TTJ 


H = 


1 f d 2 \ 

- ( + W 2 (x) + 2v + nha 3 W' J , W(x) = w 2 x 2 + w\x + w 0 . 


(16) 


For any other form of the superpotential, the nilpotent operators (pj5|) are not conserved. 
The family of supersymmetric systems (0) is reduced to the superoscillator at w 2 = 0 with 
the associated exact n-supersymmetry ||. For w 2 ^ 0, the n-supersymmetry is realized 
always in the spontaneously broken phase since in this case the supercharges (15]) have no 
zero modes (normalized eigenfunctions of zero eigenvalue). 

One can also look for the supercharges in the form of polynomial of order n in the 
oscillator-like operators A ± defined in ( 0 ); 


n —1 




(17) 


k =0 


where qu are real parameters which have to be fixed. As in the case of the supercharges 
(|15|) , the requirement of conservation of ([17]) results in the Hamiltonian ([16]) but with the 
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exponential superpotential [11]: 


W(x) = w + e ux + w-e~ ux + w 0 , 


u) 2 = 


24q 


n {n 2 — 1) ’ 


(18) 


where all the parameters w±$ are real, while the parameter u is real or pure imaginary 
depending on the sign of the real parameter q, and for the sake of simplicity we put H — 1. In 
the limit u —> 0 this superpotential is reduced to the quadratic form m via the appropriate 
rescaling of the parameters w± t0 . 

The family of n-supersymmetric systems given by the superpotential ( |T8D is tightly related 
to the so called quasi-exactly solvable problems [15|, 0, □. 0 ,0. Indeed, the both of 
the Hamiltonians constituting the supersymmetric Hamiltonian of the form (|l6j) with the 


exponential superpotential belong to the s/(2,M) scheme of one-dimensional QES systems 
||T5| , |T6[ , |T7| . Besides, the QES family given by the superpotential (|T8| ) is related to the exactly 
solvable Morse potential for some choice of the parameters 0- 

The n = 2 non-holomorphic supersymmetry corresponding to Eqs. O, (B occupies 
an especial position. Like the linear supersymmetry, it admits the anomaly-free quantum 
formulation in terms of an arbitrary superpotential. Indeed, the quantization of the super- 
symmetric system ([14]) leads to 0 


H = 


~ + v + 2hW ' a 3 + A W 

dx 2 W 2 


q + = (Q~y =; 


h l + w ) + W 2 ~^ w) 


0 + , 


A = 


h 2 

4W 2 


(2W"W - W' 2 ). 


(19) 

( 20 ) 
( 21 ) 


Looking at the quantum Hamiltonian fll9]) and supercharges (|20| ), we see that the presence 
of the quadratic in h 2 term (|2l|) in the operators H and Q + is crucial for preserving the 
supersymmetry at the quantum level. Therefore, one can say that the quantum correction 
(21) cures the problem of the quantum anomaly since without it the operators Q ± would not 
be the integrals of motion. The supercharges (^) satisfy the relation { Q + , Q~} = (H—v) 2 +c, 
and the structure of the lowest bounded states in the cases c > 0, c < 0 and c = 0 for v — 0 
is reflected in the table and on the figure (for the details see Ref. 
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Table. The structure of the lowest states for the n = 2 supersymmetry. 



c > 0 

c = 0 

c < 0 

a ) 

Completely broken phase, there are no singlet states 

+ 

+ 

+ 

b) 

One singlet state in either bosonic or fermionic sector 


+ 

+ 

c) 

Two singlet states with E — 0, one is in fermionic 

sector, another is in bosonic sector 


+ 

+ 

d) 

Two singlet states in one of two sectors 



+ 


From this structure one can see, in particular, that the quantum theory “remembers” 
its classical origin: the case c > 0 corresponding classically to the holomorphic n — 0 
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°b) 

c) 


• • 

• • 

• • 

O 

■G - 

d) 


Figure. The four types of the spectra for the n — 2 supersymmetry for 
bounded 

states. 


supersymmetry gives the systems in the completely broken phase for any superpotential 
providing the existence of bounded states. 

In conclusion of the discussion of the nonlinear supersymmetry for the ID quantum 
systems, we note that for the Erst time the close relationship between the nonlinear su¬ 
persymmetry and QES systems was observed in Ref. 0- Recently, it was demonstrated 
|32| that the so called type A TV-fold supersymmetry |33|] being a generalization of the one¬ 
dimensional holomorphic supersymmetry is, in essence, equivalent to the one-dimensional 
QES systems associated with the sl( 2, M) algebra. 


Nonlinear supersymmetry on plane in magnetic field 


The nonlinear holomorphic supersymmetry we have discussed has a universal nature due to 

m 


the algebraic construction underlying it and revealed in Ref. 


This universality allows 


us, in particular, to generalize the above analysis to the case of the two-dimensional systems. 

The classical Hamiltonian of a charged spin-1/2 particle (—e = m = 1) with gyromagnetic 
ratio g moving on a plane and subjected to a magnetic held -B(x) is given by 


H — -V 2 + gB(x)0 + r, 


( 22 ) 


where V = p + A(x), A(x) is a 2D gauge potential, B(x) = d\A 2 — d 2 A 1 . The variables 
x i} Pi , i = 1,2, and complex Grassman variables 9 ± , (6> + )* = 9~, are canonically conjugate 
with respect to the Poisson-Dirac brackets, = S t j , {9~,6 + } ir = —i. For even values 

of the gyromagnetic ratio g = 2 n, n G N, the system (|22|) is endowed with the nonlinear 
n-supersymmetry. In this case the Hamiltonian fl22| ) takes the form 

Hr, = l -Z + Z~ + l -n {Z-, Z + }* 9 + 6~ , Z ± = V 2 T iV i, (23) 


which admits the existence of the odd integrals of motion 


Q ± = 2 


(24) 
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generating the nonlinear n-superalgebra (10). The n-superalgebra does not depend on the 
explicit form of the even complex conjugate variables Z ± . Therefore, in principle, Z ± can 
be arbitrary functions of the bosonic dynamical variables of the system. 

The nilpotent quantity N = 6 + 6~ is, as in the ID case, the even integral of motion. 
When the gauge potential A(x) is a 2D vector, the system (|23|) possesses the additional even 
integral of motion L = e^XiPj. The integrals N and L generate the U( 1) rotations of the 
odd, 9 ± , and even, Z , variables, respectively. Their linear combination J = L + nN is in 
involution with the supercharges, {J, Q ± }*=0, and plays the role of the central charge of 
the classical n-superalgebra. As we shall see, at the quantum level the form of the nonlinear 
n-superalgebra ([hj) is modified generically by the appearance of the nontrivial central charge 
in the anticommutator of the supercharges. 

A spin-1/2 particle moving on a plane in a constant magnetic field represents the sim¬ 
plest case of a quantum 2D system admitting the nonlinear supersymmetry. Such a system 
corresponds to the n-supersymmetric quantum oscillator [[12]. As in the case of the one¬ 
dimensional theory, the attempt to generalize the n-supersymmetry of the system to the 
case of the magnetic field of general form faces the problem of quantum anomaly. The 
generalization is nevertheless possible for the magnetic held of special form [|l^ . 


To analyse the nonlinear n-supersymmetry for arbitrary n € N, it is convenient to intro¬ 
duce the complex oscillator-like operators 


Z = d + W{z,z), 


Z = -d + W(z,z), 


(25) 


where the complex superpotential is defined by Re ID = A 2 (x), Im W = A^x), and the 
notations z = \{x i + ix 2 ), z = \{x\ — ix 2 ), d = d z , d = d z are introduced. 

The magnetic held is defined by the relation [Z,Z] = 2 B(z,z). The n-supersymmetric 
Hamiltonian has the form 

L — 77/ — 

H n = -{Z,Z} + -[Z,Z]a 3 . (26) 

For n = 1 we reproduce the usual supersymmetric Hamiltonian. Unlike the linear supersym¬ 
metry, the nonlinear holomorphic supersymmetry exists only when the operators ([25]) obey 
the relations 


[z, [z, [Z,Z]]] =co 2 [Z,Z], 


[Z, [z, [z,z]]] =o 2 [z,z]. 


(27) 


Here u G C and u = iv*. Using Eq. (|27|), one can prove algebraically by the mathematical 
induction that for the system (E6|) the odd operators defined by the recurrent relations 


(28) 


« + 2 = i(z 2 -( 2 y 1 ) 2 ‘ j2 ) Q t. Ct = » + . Qt = 2-iZ6\ 


are the integrals of motions, i.e. they are supercharges. One can make sure [ FT | that in 
the ID case these operators generate the nonlinear supersymmetry with the polynomial 
superalgebra ([I]). 

In the representation (^5[) the conditions (|27| ) acquire the form of the differential equations 
for magnetic held: 


( d 2 — uj 2 ) B{z, z) = 0, 


(d 2 -Cj 2 ) B(z, z) = 0. 


(29) 
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The general solution to these equations is 


B(z, z) = w+e“ z+m + + toe"'-"' + ffie- 1 "'-"'), 


(30) 


where w± G R, w G C, w = w* 
the polynomial, 


On the other hand, for u = 0 the solution to Eq. (p9|) is 


B(x) = C ((xi - Xio) 2 + (x 2 - X20) 2 ) + Co, 


(31) 


with c, c 0 , (Cio, X 20 being some real constants. Though the latter solution can be obtained 
formally from (|30l) in the limit oj —> 0 by rescaling appropriately the parameters w±, w, the 
corresponding limit procedure is singular and the cases (^0|) and ([EL]) have to be treated 
separately. 

Since the conservation of the supercharges is proved algebraically, the operators Z, Z can 
have any nature (the action of Z, Z is supposed to be associative). For example, they can 
have a matrix structure. With this observation the nonlinear supersymmetry can be applied 
to the case of matrix Hamiltonians [[341, [35| , |36j| . 

Thus, the introduction of the operators Z , Z allows us to reduce the two-dimensional 
holomorphic n-supersymmetry to the pure algebraic construction. It is worth noting that in 


the literature the algebraic relations ([27]) are known as Dolan-Grady relations. The relations 
of such a form appeared for the first time in the context of integrable models |T3| . 

The essential difference of the n-supersyrnmetric 2D system (^6[) from the corresponding 
ID supersymmetric system is the appearance of the central charge 


L — — — Ti — 

J n — — ~ (cu 2 Z 2 + u 2 Z 2 ) + dBZ + dBZ - B 2 + - ddBa 3 , 


(32) 


[H n , J n ] = [Qn, Jn] =0. The anticommutator of the supercharges contains it for any n > 1. 
For example, the n = 2 nonlinear superalgebra is 


{<?2 ■ Qi } - h% + ~j 2 + Li 


(33) 


The systems (|^) with the magnetic held (^0) of the pure hyperbolic (w = 0) or pure 
trigonometric ( w± = 0) form can be reduced to the one-dimensional problems with the 


nonlinear holomorphic supersymmetry [12 


Let us turn now to the polynomial magnetic held (|3ll). One can see that this case 
reveals a nontrivial relation of the holomorphic n-supersymmetry of the 2D system to the 
non-holomorphic ID A/"-fold supersymmetry of Aoyarna et al |33 . 

In the system (]2 (d|) with the polynomial magnetic held ( |5T| ) the central charge has the 
form 


J n — (dB(z, z)Z + dB(z, z)Z — B 2 (z, z) + ^ ddB(z , z 


(34) 


It can be obtained from the operator (p2|) in the limit u —> 0 via the same rescaling of the 
parameters of the exponential magnetic held which transforms fl30f) into (^1]). The essential 
feature of this integral is its linearity in derivatives. 
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The polynomial magnetic field m is invariant under rotations about the point (x w , £ 20 )• 
Therefore, one can expect that the operator ([b|) should be related to a generator of the axial 
symmetry. To use the benefit of this symmetry, one can pass over to the polar coordinate 
system with the origin at the point (£io,^ 2 o)- Then the magnetic field is radial, B(r) = 
cr 2 + Co- The supercharges have the simple structure: Q+ = 2” 2 Z n 6 + = (Q~Y. As in the 
case u 7 ^ 0 , the anticommutator of the supercharges is a polynomial of the n-th degree in 
H n , {Qn ’ Qn } = A” + P(H n , J n ), where P(H n , J n ) denotes a polynomial of the (n — l)-th 
degree. For example, for n = 2 one has 

{ Q 2 1 Qt} = + cJ 2 . 


For the radial magnetic field it is convenient to use the gauge 

A = ^cr 4 + ^c 0 r 2 , A r = 0. (35) 

In this gauge the Hamiltonian (^) reads 

I 77 

H n = - - (d 2 r + r~ 1 d r - r~ 2 (A 2 v (r) - 2 iA v {r)d v - dfy) + -B(r)a 3 , (36) 

2 

while the central charge (|34D takes the form J n = — id^ ^ + |< 7 3 . Thus, the integral J n 
is associated with the axial symmetry of the system under consideration. The simultaneous 
eigenstates of the operators H n and J n have the structure 




e i{m-n) VXm ( r )\ 
e im *ipm{r) ) 


(37) 


Since the angular variable (p is cyclic, the 2D Hamiltonian (|36]) can be reduced to the 
ID Hamiltonian. The kinetic term of the Hamiltonian (^j) is Hermitian with respect to the 
measure dp = rdrdip. In order to obtain a one-dimensional system with the usual scalar 
product defined by the measure dp = dr , one has to perform the similarity transformation 
H n —¥ UH n U ~ l , T —> f/T with U = \fr. Since the system obtained after such a transforma¬ 
tion is originated from the two-dimensional system, one should keep in mind that the variable 
r belongs to the half-line, r G [0, 00 ). After the transformation, the reduced one-dimensional 
Hamiltonian acting on the lower (bose) component of the state (|37]) reads as 


77i 2) = - + V r4 + l ( c o - 2c ( 2n - m )) r ' + 


m 


2 r 2 


1 , 

— In 
2 v 


m)c 0 . (38) 


This Hamiltonian gives the well-known family of the quasi-exactly solvable systems 0 , 0 

T| , P^[ | . The superpartner Pin } acting on the upper (fermi) component of the state (0) can 
( 2 ') 

be obtained from Pin ; by the substitution n —> —n, m —>■ m — n. 

The reduced supercharges have the form 


a+ = 2-f ZJ + = (Q, 


-ft 


Z n = (A- A) (A-^)...A, 
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where A = + W(r) and the superpotential is 


1 1 

W (r) = -cr 3 + -c 0 r + 


m — 


The operators Tin , i — 1,2, generate the non-holomorphic type A A/"-fold supersymme¬ 


try discussed in |h|. The supersymmetry is exact for c > 0 (c < 0) and corresponding zero 
modes of the supercharge <2+ (<2“) can be found. The relation of the A/"-fold supersymmetry 
with the cubic superpotential to the family of QES system (|38|) with the sextic potential was 
also noted in Ref. 


Resume 


To conclude, let us summarize the main results of our consideration of the nonlinear super- 
symmetry. 

• Generalized statistics and supersymmetry are intimately related. 

• Linear supersymmetry at the classical level is a particular case of a classical supersym¬ 
metry characterized by the Poisson algebra being nonlinear in Hamiltonian. 

• Any classical ID supersymmetric system is symplectomorphic to the supersymmetric 
system of the canonical form characterized by the holomorphic supercharges. There 
are three different classes of the classical canonical supersymmetric systems defined by 
the behaviour of the superpotential. 

• The anomaly-free quantization of the classical ID holomorphic n-supersymmetry is 
possible for the quadratic and exponential superpotentials. 

• The nonlinear supersymmetry is closely related to the quasy-exactly solvable systems. 


The n = 2 supersymmetric Calogero-like systems (|T4| ) admit the anomaly-free quanti¬ 
zation for any superpotential; the specific quantum term ( rv./ h 2 ) “cures” the quantum 
anomaly problem. 


• The anomaly-free quantization of the classical 2D holomorphic n-supersymmetry Exes 
the form of the magnetic held to be the quadratic or exponential one. 

• Realization of the holomorphic n-supersymmetry in 2D systems leads to the appearance 
of the central charge entering nontrivially into the superalgebra. 

• The holomorphic nonlinear supersymmetry can be related to other known forms of 
nonlinear supersymmetry via the dimensional reduction procedure. 

• There is the universal algebraic foundation associated with the Dolan-Grady relations 
which underlies the holomorphic n-supersymmetry. 

The universal algebraic structure underlying the holomorphic nonlinear supersymmetry 
opens the possibility to apply the latter for investigation of the wide class of the quantum 
mechanical systems including the models described by the matrix Hamiltonians, the models 
on the non-commutative space, and integrable models 
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